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We consider a closed quantum system, initially at thermal equilibrium, driven by arbitrary exter- 
nal parameters. We derive a lower bound on the entropy production which we express in terms of 
the Bures angle between the nonequilibrium and the corresponding equilibrium state of the system. 
The Bures angle is an angle between mixed quantum states and defines a thermodynamic length 
valid arbitrarily far from equilibrium. As an illustration, we treat the case of a time-dependent 
harmonic oscillator for which we obtain analytic expressions for generic driving protocols. 
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Thermodynamics provides a generic framework to de- 
scribe properties of systems at or close to equilibrium. 
On the other hand, for systems which are far from 
equilibrium, that is beyond the linear response regime, 
no unified formalism has been developed so far. Re- 
cently, however, a number of cold atom experiments have 
been able to investigate quantum processes, which oc- 
cur far from thermal equilibrium [lHj]; they underline 
the need for general characterizations of quantum pro- 
cesses that take place beyond the range of linear response 
theory. In thermodynamics, nonequilibrium phenomena 
are associated with a non-vanishing entropy production, 
(S) = AS- (Q) /T > 0, defined as the difference between 
the change of entropy and the (mean) heat divided by 
temperature 0, 0] . The positivity of the (mean) entropy 
production is an expression of the second law of thermo- 
dynamics and follows from the Clausius inequality. The 
entropy production (E) is expected to be larger the fur- 
ther away from equilibrium a system operates. However, 
it is not possible to compute (E), nor to derive a use- 
ful, process-dependent, lower-bound for it, within equi- 
librium thermodynamics. 

For classical systems near equilibrium, such a lower- 
bound has been obtained using a geometric approach, 
and expressed in terms of the thermodynamic length 
0-0] . The latter defines a thermodynamic Riemannian 
metric which measures the distinguishability of equilib- 
rium and nonequilibrium distributions (lol | . In the linear 
regime, the entropy production is bounded from below 
by the square of the thermodynamic length (E) > £ 2 . 
The thermodynamic length plays an important role in 
finite-time thermodynamics, where it provides limits on 
the efficiency of thermal machines [TJ EH • Methods on 
how to measure £ have been discussed in Refs. [HI, H3 |. 
Interestingly, the length £ is identical to the statistical 
distance introduced by Wootters to distinguish two pure 
quantum states 15|: the angle in Hilbert space between 
two wave vectors ip\ and ip2 is given by £ (ipi,^) = 
arccos J dx \Jp\ ( x ) Vi with the two probability dis- 
tributions pi(x) = |^i(x)| 2 and P2(x) = \ip2(x)\ 2 . Re- 
cently, we have extended the notion of thermodynamic 
length to closed quantum systems driven arbitrarily far 



from equilibrium [hq . To this end, we have general- 
ized the length £ by the Bures angle C fl7l-l2l| between 
the nonequilibrium and the corresponding equilibrium 
density operators of the system. The Bures metric is 
a generalization of Wootters' metric to mixed quantum 
states and play s a major role in quantum information 
theory [22I. l23l|. Using the Bures angle, we have de- 
rived a generalized Clausius inequality with a process- 
dependent lower bound, (E) > (8/-7T 2 ) £ 2 , that is valid for 
arbitrary nonequilibrium driving beyond linear response. 
This bound, however, corresponds to the lowest order 
term of a systematic series expansion as a function of the 
Bures length. Our aim in this paper is to extend our 
previous findings and derive a sharper lower bound on 
the entropy production (E) by evaluating the contribu- 
tion of higher order terms. We then apply this result to 
the case of a quantum parametric harmonic oscillator, 
a model for a driven trapped ion [24|, [25| , for which we 
find exact analytical expressions for the angle C for ar- 
bitrary driving protocol. We furthermore compare these 
results with those obtained with the trace distance, a 
non- Riemannian quantum metric (22L [23^ . Finally, we 
derive an upper bound for the quantum entropy produc- 
tion in the Appendix. 



I. GEOMETRIC ANGLE BETWEEN MIXED 
QUANTUM STATES 

The Bures angle C is implied by the Bures metric, 
which formally quantifies the infinitesimal distance be- 
tween two mixed quantum states described by the den- 
sity operators p and p+dp as £ 2 (p+Sp, p) = tr {SpG} /2, 
where the operator G obeys the equation pG + G p = Sp 
[23j . In the orthonormal basis \i) that diagonalizcs 
p = Pi\i) an explicit expression of the Bures metric 

is given by £ 2 09+^, p) - (1/2) I (¥p\j) \V(Pi+Po)- 
In the limit of pure quantum states, the Bures metric 
reduces to Wootters' statistical distance, £ 2 (p,p + dp) = 
(1/4) J2i(dpi) 2 1 Pi [HJ. Wootters' distance is equal to the 
angle in Hilbert space between two state vectors, and is 
the only monotone, Riemannian metric (up to a constant 



2 



factor) which is invariant under all unitary transforma- 
tions [iff . It is therefore a natural metric on the space of 
pure states. The Burcs metric, on the other hand, being 
the generalization of Wootters' metric to mixed quantum 
states represents a natural, unitarily invariant Rieman- 
nian metric on the space of impure density matrices [23| . 

For any two density operators p\ and /?2, the finite 
Bures angle C is given by 



arccos 



\/F( Pl ,p 2 ) 



(1) 



where the fidelity F is defined for an arbitrary pair of 
mixed quantum states as (l9l I20I ] , 



F(pi,P2) 



tr 



Pi P2 VPl 



(2) 



The fidelity is a symmetric, non-negative and unitar- 
ily invariant function, which is equal to one only if the 
two states, p\ and p%, are identical. For pure quantum 
states, pi = \4>i)(4>i\i the fidelity reduces to their overlap, 
F(pi,Pi) = tr{pi/o 2 } = I (ipi\ip2) I 2 - It is worth empha- 
sizing that the Bures angle (|TJ) is the natural distance 
quantifying the distinguishability of two density opera- 
tors. We shall use this property in the following to quan- 
tify the distance between a noncquilibrium state and the 
equilibrium state corresponding to the same configura- 
tion of the system. With the help of this thermodynamic 
length, we will also obtain a lower bound on the nonequi- 
librium quantum entropy production. 



II. THERMODYNAMIC LENGTH AND 
GENERALIZED CLAUSIUS INEQUALITY 

We consider a quantum system whose Hamiltonian 
H = H t is varied during a finite time interval r. We 
assume that the system is initially let to equilibrate with 
a thermal reservoir at inverse temperature /3 = be- 
fore an external control parameter is modified. We fur- 
ther assume that the system is quasi-isolated during the 
finite driving time r, so that relaxation is negligible and 
the dynamics is unitary to an excellent approximation. 
This corresponds to a realistic experimental situation. 
For an infinitely large driving time, much larger than 
the relaxation time induced by the weak coupling to the 
reservoir, the transformation is quasistatic and the sys- 
tem remains in an equilibrium state at all times. During 
such a slow, equilibrium transformation, the change in 
free energy, AF = AE — TAS, is equal to the average 
work (W) done on the quantum system, (W) = AF. 
Here AE = (H T ) — (Hq) is the (internal) energy differ- 
ence. For a fast, noncquilibrium transformation, work 
is larger than the free energy difference. Using the first 
law, AE = (W) + (Q) , we can rewrite the nonequilibrium 
entropy production (E) = AS — (Q) /T as, 



The noncquilibrium entropy production (E) is thus pro- 
portional to the difference between the noncquilibrium 
and the equilibrium work done on the system. Being 
a mechanical quantity, it is worth noticing that work is 
always defined, even for far from equilibrium processes. 

Let us denote the density operator of the system at 
time t by pt- The initial equilibrium density operator is 
then p = exp(~f3H )/Z , where Z = tr {exp(— (3H )} 
is the initial partition function, with similar expressions 
for the equilibrium density operator p° q and the parti- 
tion function Z T at the final time r. To obtain a mi- 
croscopic expression for the entropy production, we use 
AE = ti{p T H T } — tr{p H n } and note that —f3H n ^ T = 
In po.r + In Zq_ t . Combined with the expression — /3AF = 
— \n(Z T /Zo) for the free energy difference, we find 

(E) = S (p T \\p^) = tr{p T ]np T - p T lnp?} , (4) 

where S(p T \ \p% q ) is the quantum relative entropy [26l[27|. 
Equation ^ is an exact expression for the nonequi- 
librium entropy production for closed quantum systems 
driven by an external parameter, and a quantum gener- 
alization of the classical results presented in Refs. (28l.[29j 
(see also Ref. H3|). We note, however, that the relative 
entropy is not a true metric, as it is not symmetric and 
does not satisfy the triangle inequality; it can therefore 
not be used as a proper quantum distance [3l[ . We next 
derive a lower bound for the quantum entropy production 
which we express in terms of the Bures angle (JTJ) . 

Inequalities are important tools of classical and quan- 
tum information theory, as they allow to express 'impos- 
sibilities', that is things that cannot happen [3l|. An 
elementary example is Klein's inequality, S(p\\\p2) > 0, 
which expresses the non-negativity of the quantum rel- 
ative entropy [22|. Combined with Eq. it immedi- 
ately leads to the usual Clausius inequality. A general- 
ized Clausius inequality can be derived by noting that the 
quantum relative entropy satisfies (Ref. [32| , Theorem 4) , 



S(pi\\p2) > s 



d(pi,P2) 

1.1 t>2,2 



d(e 



(5) 



if d(pi, P2) is an unitarily invariant norm. Further, 
\i)(j\ is the matrix with i,j elements equal to 1 and all 
other elements 0. The lower bound © has been derived 
with the help of optimization theory and is therefore as 
sharp as possible. The function s(x) is explicitly given 
by the expression [32| , 



i(x) = min I (1 — r + x) log (1-1 ) 

x<r<l I \ 1 — r J 

+ (r-a;)log(l-~)}. 



(6) 



(Z) = 0((W) - AF). 



(3) 



The first five nonzero terms in a series expansion around 
the origin x = read, 

9 4 4 32 fi 992 o 6656 in 12s 

s(x) = 2x 2 + -x 4 + x 6 + x 8 + x 10 + O(x 12 ) . 

w 9 135 5103 32805 v ; 

(7) 
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Applying inequality © to the unitarily invariant Burcs 
angle C, we obtain a process-dependent lower bound on 
the noncquilibrium entropy production. Taking into ac- 



count that C(e 



*) = it/2, since the two matrices e 



and e 2,2 are orthogonal [F(e 1,1 ,e 2,2 ) = 0], we find, 



C{p 7 



p e r q : 



> 



C 2 (p T ,p T q ). (8) 



The first order term in the expansion ((?]) yields the gen- 
eralized Clausius inequality, (£) > (8/ir 2 )£ 2 , derived in 
Ref. [l6|. Since the terms in the expansion ([7]) are posi- 
tive, an increasingly sharper lower bound can be obtained 
by taking more terms into account [33| ■ An illustration 
for the case of a quantum harmonic oscillator with time- 
dependent frequency, to be discussed in detail in the next 
Section, is shown in Fig. 1. 

Equation ([S]) indicates that the nonequilibrium entropy 
production (E) is bounded from below by a function of 
the geometric distance between the actual density oper- 
ator p T of the system at the end of the driving and the 
corresponding equilibrium operator p^ q , as measured by 
the Bures angle. Thus the Bures angle provides a natu- 
ral scale to compare (£) with, and quantifies in a precise 
manner the notion that the entropy production is larger 
when a system is driven farther away from equilibrium. 

In the classical limit, where nonequilibrium and equi- 
librium states are diagonal in the energy basis, the Burcs 
angle reduces to Wootters' statistical distance. As a re- 
sult, Eq. ([5]) yields a lower bound to the classical noncqui- 
librium entropy production that is valid for any noncqui- 
librium driving beyond the linear response regime, 



1 



(E) cl > S (-f(p T ,i^) 



> 



i 2 {p T ,p T q ). (9) 



Moreover, when nonequilibrium and equilibrium states 
are infinitesimally close, Eq. ([9]) takes the form (E) cl > 
(2/ir 2 )d£ 2 , which has been obtained in Rcfs. It 
is worth emphasizing that the latter was derived by ex- 
panding the entropy around equilibrium to second order; 
it is therefore only valid in the linear response regime. 



III. PARAMETRIC HARMONIC OSCILLATOR 

Let us now apply the generalized Clausius inequality 
© to the case of a time-dependent harmonic oscilla- 
tor. The latter provides an important physical model for 
many quantum systems, for example ultracold trapped 
ions [Ij, [2{| , and is furthermore analytically solvable. We 
will, in particular, evaluate the Bures angle ([1]) for a non- 
trivial quantum time evolution. Explicit expressions for 
C are in general only known for low dimensional systems 
[34U3H ] . The difficulty arises from the operational square 
roots in the definition of the quantum fidelity @. For 
Gaussian states, however, the expression for the fidelity 
simplifies and can be written in closed form [37j . 



The Hamiltonian of the parametric quantum harmonic 
oscillator is of the usual form (M denotes the mass), 



H t = 



P 

2M 



M 



(10) 



We assume that the time-dependent frequency u) t starts 
with initial value ojq at t = and ends with final value oj\ 
at t = t. Due to the quadratic form of the Hamiltonian 
(fTUjl . the wavefunction of the oscillator is Gaussian for 
any driving protocol uj±. By introducing the Gaussian 
wave function ansatz 13811 . 



ipt{%) = exp 



2h 



(a t x 



2bfX + Ct] 



(11) 



the Schrodinger equation for the quantum oscillator can 
be reduced to a system of three coupled differential equa- 
tions for the time-dependent coefficients at, bt and Ct, 



M 

d t b t 
d t c t 



ih 



± 

M 2 

^b 
M 

(it 



1 



M M 



(12a) 
(12b) 

(12c) 



The nonlinear equation (|12a[) is of the Riccati type. It 
can be mapped onto the equation of motion of a clas- 
sical, force free, time-dependent harmonic oscillator via 
the transformation a t = MX t /X t . The resulting equa- 
tion reads X t + uj 2 X t = 0. Equations (|12b[) - (|12cp can 
be solved once the solution of Eq. ()12ap has been deter- 
mined. With the solutions of the three equations (|12a|) - 
(|12cp known, the Gaussian wave function ipt(x) (fTT|) is 
fully characterized by the time-dependence of the angu- 
lar frequency Ut- It can be shown that the dynamics is 
completely determined by the function Q* introduced by 
Husimi 



2wq w i 



■ x i + x: 



Yr + Y: 



(13) 

where X t and Y t are the solutions of the force free clas- 
sical oscillator equation satisfying the boundary condi- 
tions Xq = 0, Xq = 1 and Yq = 1, Yq = 0. The function 
Q* > 1 is a measure of the adiabaticy of the process: it is 
equal to one for adiabatic transformations and increases 
with the degree of nonadiabaticty. In particular, the final 
mean energy of the quantum oscillator is given by [40l | , 



coth(/3fiwo/2), 



(14) 



and thus linearly increases with Q*. 

To evaluate the Bures angle (|T|) for the parametric har- 
monic oscillator, the quantum fidelity (|2|) has to be writ- 
ten in closed form. For Gaussian states such an explicit 
form is known: for two arbitrary (non-displaced) Gaus- 
sian density operators, p\ and P2, the fidelity reads |37j . 



F(pi,P2) 



(15) 
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The two parameters A = det {A\ + A2) and 6 = 
(det (Ai) — 1) (det (A 2 ) — 1) are completely determined 
by the covariance matrices A4 (matrices of the variances 
of position and momentum) of the quantum oscillator, 



(16) 

(17a) 
(17b) 
(17c) 



The matrix elements a 1 are explicitly given by 

<4 = 2 ((*?) - <^> 2 ) 

a PP = J£ ((Pi) - (Pif) 
2 /l 

(XiPi +PiXi) - (Xi) (pi) 



Xp 



h\2 



To evaluate the terms appearing in the Clausius inequal- 
ity (HJ), we make use of the explicit expressions of the 
initial, poi and final density operators, p T and p° q , of 
the oscillator in coordinate representation, as given in 
Appendix C. In particular, the final equilibrium density 
operator p° q has the same form as Eq. (|C1[) , replacing loq 
by uj\ . Accordingly, the corresponding equilibrium mean 



and variances are 



(xX 



= (p) C T q = ( X P + P x Tt = 0, and 



2Mwi 



/ 9x cq fajJ\M 



coth (fihui/2) 
coth 08fiwi/2) 



(18a) 

(18b) 
I 



On the other hand, for the final noncquilibrium state p T , 
we have (x) T = (p) r = 0, and 



(P 2 ) T 



2Mw 

hM 

2uj 



Y* + Lj£Xt) coth (/3fiwo/2)(19a) 



*+u%X*) coth(/3fi^ /2). (19b) 



r 2 



The cross correlation function can be evaluated by ex- 
ploiting the fact that (xp + px) T = M d t (x 2 ) T , and reads 



(xp + px) r = — Y T Y T + cu 2 X T X T coth (/3huj /2) . 
wo V / 

(20) 

The analytic expression of the quantum fidelity function 
F(p T ,p° q ) between noncquilibrium and equilibrium os- 
cillator states at the end of the driving can be finally 
obtained by evaluating the determinants A and S in 
Eq. ([T5J using Eqs. ([Tga |) -([2"0" |) . with the help of the re- 
lation X t Y t — X t Y t — 1 [H| and the definition of the 
function Q* given in Eq. (Til?)) . We find, 



y/ci 2 (f3e /2) + ct 2 (/3 £l /2) + 2Q* ct(/3e /2)ct(/3 ei /2) 

I 



c 2 (/3 eo /2)c 2 (^ ei /2) - c(/3e /2)c(/3ei/2) ' 



(21) 



with the notation c(.) = csch(.) and ct(.) = coth(.), and 
the energies £j = fkJi. The Burcs angle C then directly 
follows from Eq. (TTJ) , and the lower bound to the noncqui- 
librium entropy production (E) can be determined, to 
any order, with the help of the expansion ([7]). 

To get more physical insight, let us evaluate the limit- 
ing expressions of the fidelity (f2"Tj) in the low-temperature 
(quantum) and high-temperature (classical) regimes. An 
expansion of the hyperbolic cosine and cotangent func- 
tions in the zero-temperature limit, Hj3 — > 00, leads to, 



F(Pr, P C T q ) 



l + Q* 



(22) 



In the adiabatic limit Q* — > 1, the fidelity thus tends 
to one, that is the Bures angle approaches zero, indi- 
cating that the system ends in an equilibrium state, as 
expected. For strongly nonadiabatic processes, Q* >■ 1, 
on the other hand, the fidelity tends to zero as 1/yfQ*. 
Here the Bures angle tends to tt/2, showing that p T and 
p^ q are maximally distinguishable (orthogonal). 



Equation (|2"2"j) can also be derived directly by noting 
that in the zero-temperature limit, the harmonic oscilla- 
tor is initially in a pure state | Oo) - The initial, equilib- 
rium density operator is, hence, po|r=o = I Oo) (Oo | and 
analogously for p° q | T Q . Since these states are pure, 
the fidelity simplifies to their overlap, and we have, 
F (Pr, Pr q ) |t=o = ti{PrP c T q } = Po fi i wh erc p^ is the 
probability for the system to start and end in the corre- 
sponding ground state. The latter is given by the expres- 
sion H, 



Po.o 



1 



(23) 



and we thus recover Eq. 

In the classical limit, H/3 — > 0, by repeating the same 
analysis, the fidelity (f2Tj) simplifies to, 



F(p T ,p?) h n° 



4 LVqUJi 



ujq + 2 Q*uj oji + u>{ 



(24) 
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For an adiabatic frequency change, Q* —> 1, the fidelity 
reduces to F (p T , p% q ) ~ 4oj wi/(a;o + wi) 2 . Therefore, 
as noticed in Ref. [39| (see also Ref. [4l[), a unitary pro- 
cess can only be quasistatic in the thermodynamic sense, 
F (p T , p° q ) ~ 1, if \ui — wo |/^o -C 1. In addition, we note 
that in the classical limit the fidelity vanishes for large 
values of Q* as 1/Q*, that is much faster than in the 
low-temperature regime. The density operators p T and 
p eq thus become orthogonal (C = 7r/2) much faster as a 
function of the degree of nonadiabaticity. 

For the parametric quantum oscillator, the nonequilib- 
rium entropy production §5§ can be determined exactly, 
allowing to test the generalized Clausius inequality ©. 
It is given by [39|, 

(E) =| (Q*ftwi - fiw )coth(^fto;o/2) 

ln / sinh (gfafr/2) \ ( 25) 
VsmH8^V2) J ' 

where we used (Ho) = (huJo/2) coth(/3fiwo/2) and 
Eq. (I14p . Figure 1 shows the nonequilibrium entropy pro- 
duction (E) as a function of the measure of adiabaticity 
Q* , together with the lower bound obtained with the first 
term in the expansion ([7]) and the exact function s(x) 
© (the latter is indistinguishable from the expression 
including the first five nonzero terms of the expansion). 
We see that the first term in the expansion provides a 
good lower bound in many cases. 




FIG. 1: (Color online) Irreversible entropy production (E) 
(|25fl (blue, solid line) together with the lower bound Q cor- 
responding to the lowest order term of the expansion ©(red, 
dashed line), or the exact function s(x) (J6]) (green, dotted 
line), as a function of the adiabaticity parameter Q* (|13|) . 
Parameters are h = 1, ft = 1.2, cjq = 0.9, uii = 0.5. 



IV. LOWER BOUND BASED ON THE TRACE 
DISTANCE 

As discussed in Section I, the Bures angle C, being 
the extension of Woottcrs' statistical distance to mixed 



states, possesses a simple interpretation as the geometric 
angle between two density operators. However, Eq. §5§ 
shows that the nonequilibrium entropy production (E) 
is bounded by many unitarily invariant distances, albeit 
with possibly less natural physical interpretation. To elu- 
cidate this point, we discuss the concrete case of the trace 
norm, which has been reported to yield the largest lower 
bound on the relative entropy [32j (a further comparison 
of C with the Bures distance T> is given in Appendix A). 
The trace distance between two density operators, p\ and 
p 2 , is defined as (HHI 

T(pi, P2 ) = |tr{| Pl - pa)} = itr{V(pi-p 2 ) 2 } (26) 

Contrary to the Bures angle (or the Bures distance), it 
is not a Riemannian distance — however, both are mono- 
tone. The trace distance between nonequilibrium and 
equilibrium states of the parametric quantum oscillator 
(|10[) can be evaluated for arbitrary driving with the help 
of the explicit expressions of p T and p oq given in Ap- 
pendix C: we have T(p T ,p c -?) = (1/2) J2i Mi where Xi 
are the eigenvalues of p T — p° q . Unlike for the Bures angle, 
it does not seem to be possible to express T as a function 
of the adiabaticity parameter Q* alone (the density op- 
erators p T and p° q depend on the two functions X t and 
Y t , and not on Q* directly). To circumvent this problem, 
we have numerically evaluated the trace distance for the 
case of a sudden switching of the frequencies for which 
Q* = (wg + uj : 1) /(2ui qui 1). Figure [2] shows the correspond- 
ing entropy production (E) (|2"5]) and the lower bound © 
for the Bures angle, s(2£/n), and for the trace distance, 
s(T), as a function of Q* for fixed /3 and ujq. Contrary 
to Fig. [I] where both loq and lu\ arc fixed, for the sudden 
frequency switch Q*(lo\) is a function of u)\. 
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FIG. 2: (Color online) Irreversible entropy production (E) 
(|25[) (blue, solid line) together with the lower bound (JSj) eval- 
uated for the Bures angle C {U (green, dotted line), and the 
trace distance T (126[l (purple, dotdashed line), as a func- 
tion of the adiabaticity parameter (|13[1 for the sudden switch 
Q*(u)i) — (ujq + ujf )/(2cjou;i). Parameters are H = 1, f3 = 4.8, 
coo = 0.9. 

We observe that the lower bound based on the trace 
distance is sharper than the one obtained using the Bures 
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angle. However, the two bounds appear largely equiva- 
lent, reflec ting the fact that C and T are closely related 
(see e.g. [22t l23l|). We stress that the trace distance 
lacks the simple interpretation of the Bures angle as the 
angle between density operators. Moreover, in the classi- 
cal limit the bound based on the trace distance does not 
reduce to the known bound on the entropy production 
derived in linear response theory ■ 



V. CONCLUSIONS 

The Bures angle between the nonequilibrium and the 
corresponding equilibrium state of a driven closed quan- 
tum system defines a thermodynamic length that is valid 
arbitrarily far from equilibrium. The latter can be used 
to characterize the departure from equilibrium for generic 
driving. We derived a lower bound on the nonequilibrium 
entropy production, which we expressed as a function of 
the Bures angle, by using a sharp lower bound on the 
quantum relative entropy. In such a way, we obtained a 
generalized Clausius inequality, with a process-dependent 
lower bound, that holds beyond the range of linear re- 
sponse theory. As an illustration, we treated the case 
of a time-dependent harmonic oscillator for which we de- 
rived analytic expressions for the Bures angle. We further 
compared the lower bound obtained with the Bures angle 
with the one based on the trace distance. While the trace 
distance offers a slightly sharper bound, the two appear 
to be largely equivalent. 
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Appendix A: Lower bound based on the Bures 
distance 



To evaluate the changes induced by the choice of an- 
other unitarily invariant distance on the lower bound © , 
we present in this Appendix an alternative, constructive 
derivation of the lowest order estimation of the nonequi- 
librium entropy production (E). Let us begin by intro- 
ducing the Hellinger distance [42| . 

H 2 (pi,p 2 ) = [ dx (yfcjx)- VMx)) 2 , (Al) 

for two (classical) probability distributions pi(x) and 
p 2 (x). The Hellinger distance is another measure of the 
distinguishability of two probability distributions. It is a 
true distance which fulfills symmetry, non-negativity and 



the triangle inequality. Expression (|A1[) can be rewrit- 
ten in terms of the classical fidelity function, f(pi,P2) = 
/ dxy/pi{x)p 2 (x), to yield 



H( Pl ,p 2 ) = V2-2/(p 1)P2 ). (A2) 
By using the inequality, y/y — 1 > 1/2 In (y), we have 



/ P*W _ i > J {lnp2{x) _ lnpi(x)) . (A3) 

V pi i x ) 2 

Averaging Eq. (|A3j) over the probability distribution 
Pi(x) results in, 



2 1 



IP2 (x) 



< (\n Pl (x)-\np 2 {x)) pi , (A4) 



V pi ( x ) / 

' 'Ply 

from which we deduce the inequality, 

D(j>i\\p 2 ) > 2-2/(pi,p 2 ) = H 2 (pi 1 p 2 )- ( A5) 

Here D(pi\\p 2 ) denotes the classical Kullback-Leibler di- 
vergence between p\ and p 2 . The classical result (|A5|) 
can be extended to quantum states by considering the 
quantum version of the Hellinger distance which is the 
Bures distance between density operators p\ and p 2 [20j , 



V 2 (p 1 ,p 2 ) = 2(l-y/F(p 1 ,p 2 ) 



(A6) 



Note the difference in the definitions of the classical and 
quantum fidelity. By combining Eqs. (|A5j) and (|A6[I . 
we then find the generalized Clausius inequality, 



(£> >V 2 {p Tl p 



cq % 



(A7) 



The above lower bound on the entropy production cor- 
responds to the lowest order term in the expansion ([7]) 
when the Bures distances is chosen instead of the Bures 
angle in Eq. ([5]), since V(e 1 ' 1 , e 2,2 ) = y/2. Figure [3] shows 
the Bures angle C and the Bures distance I? as a function 
of the quantum fidelity F. We observe that T> > C so 
that the Bures distance offers a (slightly) sharper bound 
to the entropy production than the Bures angle, to lowest 
order. However, the distance T> bears the disadvantage 
that the intuitive, physical interpretation as an angle be- 
tween mixed states is lacking. 



Appendix B: Analytic upper bound for the quantum 
relative entropy 

Due to the importance of the relative entropy in 
physics and mathematics, and the complexity to evalu- 
ate it, accurate approximations and bounds are essential. 
While a lower bound has been obtained for unitarily in- 
variant norms in the form of Eq. ([5]) [Hf , upper bounds 
are much more difficult to find. Recently, a general up- 
per bound was proposed in terms of the eigenvalues of 



7 




FIG. 3: (Color online) Lower bounds for the nonequilibrium 
entropy production based on the Bure angle, (8/n 2 )C 2 (JSJ) 
(blue solid), and on the Bures distance, T> 2 (|A7|) (red dashed), 
as a function of the fidelity F. 



the density operators [44| . In the present thermodynamic 
context, we may, however, derive a simpler upper bound. 
We start with the inequality [45j . 



tr^lnpx -pilnp 2 } < -tr {p\+» p^ v - p x } , (Bl) 

which is true for all positive definite operators p\ , p-2 
and v > 0. We shall here concentrate on the final 
nonequilibrium and equilibrium density operators, p T 
and /5® q . By choosing v = 1 and using the normaliza- 
tion condition tr {p} = 1, we obtain the upper bound, 

(£) < trj/^ (p^r 1 } - 1. In order to further simplify 

the bound and derive an expression which does not de- 
pend on the off diagonal matrix elements in energy repre- 
sentation of the density operators, we use the inequality 

m, 



|trW 2 }| < 



\ " 12 

cr„ a,. 



(B2) 



which holds for any complex n x n matrices p\ and p2 
with descending singular values, a\ > ... > a\ and 
a\ > ... > a 2 . The singular values of an operator A 



acting on a Hilbert space are defined as the eigenvalues 
of the operator v A^ A. If p\ and pi are density opera- 
tors acting on the same Hilbert space Eq. (|B2jl remains 
true for arbitrary dimensions, and the singular values are 
identically given by the eigenvalues J47[- As a result, we 
obtain the upper bound for the entropy production (E), 



< 



E 



(pIY 



(B3) 



Appendix C: Explicit expressions of the density 
operators 



The evaluation of the covariance matrix (|16|) requires 
the expressions of the density operators pg q , and p T 
in coordinate representation. We collect them in this Ap- 
pendix for completeness. The initial equilibrium density 
operator pg q is given by (49j . 



AIojq 



tanh (f3huJo/2)x 



exp 



Mu>, 



2/) 



- coth (fihujo) (x 2 + y 2 — 2sech(/3fta;o)a;y) 



(CI) 



The final equilibrium density operator has the same 
form as Eq. (|C1[) with the replacement uiq by uj\. On 
the other hand, the final nonequilibrium operator p T can 
be derived from Eq. (|C1|) by noting that p T (x,x') = 
f dy f dy' U T (x,y) p c q (y,y')U;(y',x'). The propagator 
[/ t (.t, .t ) can be determined from the wave function (fTTj) 
with i/Jt{x) = J dxoU T (x,xo)tpt (xo), and reads [38j ] - 



*7r = 



A/ 



exp 



2nihX T V 2fi.A r 



Y t Xq 



(C2) 

The functions X r and Y r are solutions of the force free 
harmonic oscillator satisfying the boundary conditions 
Xq = 0, Xq = 1 and 1q = 1, lo = 0. A direct evaluation 
of the Gaussian integrals leads to the expression, 



Pr(x,y) 



'Mcj tanh(/3Sw /2) 



exp 



Mojq 1 
2h Y?+lo 2 X* 



"I* 2 



coth(/3^ ) {x 2 + y 2 -2sech(f3hiuo)xy) + i (x 2 - y 2 ) (uj 2 X t X t + Y T 



Y T 



(C3) 
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